Noise resilience and entanglement evolution in two non-equivalent classes of quantum 

algorithms 
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The speed-up provided by quantum algorithms with respect to their classical counterparts is 
at the origin of scientific interest in quantum computation. However, the fundamental reasons 
for such a speed-up are not yet completely understood and deserve further attention. In this 
context, the classical simulation of quantum algorithms is a useful tool that can help us in gaining 
insight. Starting from the study of general conditions for classical simulation, we highlight several 
important differences between two non-equivalent classes of quantum algorithms. We investigate 
their performance under realistic conditions by quantitatively studying their resilience with respect 
to static noise. This latter refers to errors affecting the inital preparation of the register used to 
run an algorithm. We also compare the evolution of the entanglement involved in the different 
computational processes. 

PACS numbers: 03.67.-a, 03.67.Lx, 03.67.Mn 



I. INTRODUCTION 



One of the major reasons for investigating quantum 
computation is the possibility for a quantum processor 
to outperform any analogous classical device [l[ . The de- 
sign of genuine quantum algorithms and the study of the 
reasons for their speed-up with respect to classical coun- 
terparts have been the center of considerable interest. It 
is now generally accepted that, for pure states, quantum 
correlations spread over a sufficiently large number of el- 
ements of a register [31 and quickly growing with the size 
of the register itself [3(] are a necessary requirement for 
the speed-up. However, for registers prepared in mixed 
states, the requirements are still largely unknown. It has 
been conjectured that, for mixed states, the criteria men- 
tioned above are not sufficient. 

The study of classical simulation of quantum algo- 
rithms can help us in understanding the role of inherently 
quantum phenomena in computational problems. Con- 
siderable effort has been made in this direction, with pro- 
posals for classical simulations designed in experimental 
setups ranging from nuclear magnetic resonance \4 , to 
cavity-quantum electrodynamics and linear optics [5|. It 
has been pointed out by Meyer that it is possible to clas- 
sically simulate quantum algorithms that rely on the use 
of balanced linear superpositions of the computational 
states of a register Q. Such an initial state can then 
be reinterpreted as the state of a multilevel particle by 
neglecting the multipartite nature of the register. This 
allows one to reinterpret quantum entanglement in terms 
of simple coherences. In this paper, we use the quan- 
tum average algorithm Q as the representative of a class 
of quantum protocols (from now on indicated as non- 
polylocal) which is non-equivalent to the class identified 
by Meyer (labeled as polylocal). The latter is represented, 
in our study, by the quantum search algorithm 0. 

A clear difference between the two classes is the "na- 
ture" of the initial state of the register. Polylocal algo- 
rithms use initially separable states and generate entan- 



glement during their performance [9]. Differently, non- 
polylocal protocols exploit entangled initial resources. As 
a result, in this second class of algorithms, the state of 
the register cannot be put in correspondence with an un- 
biased state of a multilevel system. This prevents the 
use of general arguments a la Meyer. Moreover, the en- 
tangled resource itself is "consumed" during the process- 
ing of non-polylocal algorithms. It is also interesting to 
notice a close analogy between these protocols and the 
measurement-based model for computation [Tfj| . 



The different use of entanglement in the two classes of 
problems makes their quantitative comparison difficult. 
In our study, we use the influence of noise (and thus the 
introduction of classical correlations in the algorithms) as 
a useful tool for the investigation into differences between 
the representatives of polylocal and non-polylocal algo- 
rithms. For the quantum average algorithm, although a 
fragile GHZ-like state [ll| is used, we find a consider- 
able resilience to static noise. Moreover, the GHZ-like 
nature of the resource remains unchanged, leading us to 
conclude that this specific form of entanglement plays a 
crucial role in the performance of the algorithm. 



This paper is organized as follows. In Sec. |TT] we pro- 
vide a brief explanation of how the representatives of the 
two different classes work and briefly comment on their 
classical simulation. In Sec. lIIII we investigate the perfor- 
mances of these protocols in the presence of static noise, 
focusing our attention on the quantum average algorithm 
for which, to the best of our knowledge, noise resilience 
has never been studied. The evolution of entanglement 
is studied in Sec. IIV1 where a clear picture of the salient 
properties of the quantum average algorithm is provided. 
Finally, in Sec. [V] we summarize our results. 



II. DYNAMICS OF THE REPRESENTATIVE 
ALGORITHMS 

A. Polylocal class: Quantum search algorithm 

The quantum search algorithm [8| is designed to find a 
searched item in a randomly ordered database of length 
L in an 0(\/L) time [Hj]. If we want to carry out the 
same search using a classical algorithm, on average 
steps are required, as we need to analyze the items one 
by one until the searched one is found. In order to de- 
scribe the algorithm, we assume that each item of the 
database is labeled by a binary number between and 
L — 1 with L = 2™ and n an integer. Thus, the task 
becomes the identification of the number labeling the 
searched item. Using an n-qubit system, each state of 
the computational basis corresponds to a binary num- 
ber in the set {Q,...,L — 1}. For example, the state 
|0 • • • 0101) corresponds to 101. The algorithm con- 
sists of an alternating sequence of operators, as a re- 
sult of which the target state can be found in O(VZ) 
queries. The initial state of the register is the superposi- 
tion 1 6) = (l/\/Z) li) with \i) one of the states of the 
computational basis and i = 0, L— 1. This state can be 
obtained by applying H® n to |0) with H the Hadamard 

transform H = (\ ^~^\ written in the single-qubit 

basis {|0),|1)}. The sequence of operators, or Grover 
iterate, is made out of the following four steps: 

• Query of the oracle [l3|: The oracle performs a 
phase-flip on the searched state. This is effectively 
performed through the operator 1—2 \ s)(s\ with |s) 
the searched state and I the L x L identity matrix. 

• Application of H® n . 

• Application of 2 |0) (0| - 1 

• Application of H® n . 

The query of the oracle can be seen as being at the heart 
of the algorithm. For a generic input state, the oracle 
marks the searched state flipping the sign of its ampli- 
tude. Classically, this corresponds to a checking function 
which produces a different output depending on whether 
the input is the searched one or not. The last three 
operations can be seen as the application of the oper- 
ator 2 1 0) (0 1 — I, which carries out an inversion about 
the mean. By iterating this sequence, the state of the 
n-qubit system oscillates between the equally weighted 
superposition 1 5) and the searched state. The first max- 
imum of this oscillation is obtained after R iterations of 
the sequence, where R is the closest integer to the real 
number 
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FIG. 1: Logical circuit for the quantum average algorithm [?|. 
The input is the GHZ state in Eq. (J5]). Single-qubit rotations 
are indicated as R(ctj) with otj = jfg and the Hadamard gate 
H is also shown. The qubits are measured in the a x eigenbasis 
(the H gate and the (Jz-basis measurement can be seen as a 
a-x-basis measurement). The dashed line represents classical 
information. 



For large L, X (and therefore R) becomes (D(y/~L), so 
that the number of steps required to find the state \s) 
with almost certainty scales as \[L. 



B. Non-polylocal class: Quantum average 
algorithm 

Suppose we have ./V values Uj S [—1, 1]. A well-known 
computer-science problem is to find the order of magni- 
tude of the average /i, defined by /i = YljLi v r Clas- 
sically, if we pick up m random samples from the N- 
value set, the average evaluated out of them will be dis- 
tributed according to a Gaussian centered at the actual 
average and with standard deviation 0(^=), as obtained 
in virtue of the central limit theorem. This means that, 
with high probability, the estimated average lies within 
0(— L) of the true average. Suppose that, for a specific 
problem to be solved, /i must be known with a precision 
at least equal to e. This means that 0(-i=) < e and thus 

m > samples are needed to estimate the average 

with a precision of e. If we want to know the order of 
magnitude of the average, we need ^(^r) samples. A 
speed-up in finding the solution to this problem can be 
achived by using the following quantum algorithm, which 
is able to estimate the ratio for a fixed 9 > 0, in a 
number of steps independent of u and 6 and depending 
only on the precision we want to have with respect to the 
estimate of this ratio. 

Consider a register of N qubits prepared in the gener- 
alized GHZ state 
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(see Fig. ??) The underlying assumption is that the val- 
ues of the set {i/j} are distributed to the N stations of 
a network. Each agent has knowledge of just the value 
that has been attributed to him. We shift the phase 
of the j-th qubit (j = 1,...,N) by b y a PPly in g the 



2 arccos 



L-l 
L 



|0) (0| + |1> (1|. Af- 
ter these single-qubit operations, the state of the register 



(1) conditional operator Rj( 



3 V N9 > 



3 



becomes 

= ^(e** |0>! |0> 2 ---|0> iV + |l) 2 -|l)iv)- 

By means of cr^-basis measurements on each qubit but 
the first one, we end up with |?/' ± )i = "^f ( e * T |0) 1 ± |1) 1 ) 
with the plus {minus) sign if the number of f's mea- 
sured on the other qubits is even (odd) and a r the r- 
Pauli matrix (r = x, y, z). If the number of l's is odd, 
we can obtain simply by shifting the phase of the 

first qubit by tt. The information about (i is now car- 
ried by the first qubit and we can estimate it by an 
interference-type procedure where we apply H to obtain 

^|±i 10)! + £%i |l) r The probability of measuring 
|0) or |1) out of this state are now respectively cos 2 (^) 
and sin 2 (^). Therefore, by repeating this procedure a 
times, in virtue of the central limit theorem, we can esti- 
mate ^ with a precision 0(^=). A close analogy with 
measurement-based computation can be seen here, as the 
additional phase-shift conditioned on the outcomes of the 
previous N — 1 measurements represents the byproduct 
operator of the specific scheme at hand [ToL [l4| . 

If we need to estimate the order of magnitude of /i, 
we can start the algorithm by taking a large value of 9 
(say, for example, 0.5) and evaluate the ratio |/i|/6>. If 
this ratio is found to be 0(1), we have the correct or- 
der of magnitude of /i. Differently, if the ratio is much 
smaller than one, we divide 6 by a fixed value (say 2) 
and estimate again \^\/9 until we find it to be 0(1) 0- 
In this way, we need O(log 2 /i) applications of the quan- 
tum algorithm to solve the problem (each application of 
the algorithm requires 0(a) steps, but if we fix the pre- 
cision we want in the estimate, a is fixed too). On the 
other hand, we have seen that any classical algorithm can 
solve this problem in steps, therefore a speed-up 

is provided by the quantum algorithm. This protocol is 
particularly useful in a scenario of distributed quantum 
computation, where a processor is made out of a network 
of local nodes interconnected by classical and quantum 
channels @, [H| ■ Indeed, the qubits can be at remote lo- 
cations and can operate independently with the partial 
knowledge of just the value of the belonging to node 
j. The only requirement for the algorithm to work, once 
the initial resource is provided, is the transmission of one 
bit of classical information (the result of each measure- 
ment) to the first qubit location. 



C. Classical simulations of quantum algorithms 

The classical wave optics analogy of quantum infor- 
mation processing is based on the fact that the state of 
a quantum system evolves according to a wave equation 
and satisfies the superposition principle [l[. Many clas- 
sical simulations of quantum algorithms have been pro- 
posed in recent years. They require a number of classical 



resources scaling exponentially with the number of qubits 
being simulated [Ty] . For example, a few classical opti- 
cal simulations represent the Hilbert space of n qubits by 
considering the propagation of a classical electromagnetic 
wave. Splitting the cross section of this wave in 2™ dif- 
ferent spatial zones allows one to associate the amplitude 
of the electromagnetic wave in each zone with the ampli- 
tude of a state of the computational basis of the quantum 
system to be simulated [5|, [13] • Another proposal put 
forward is to represent n qubits by a single photon in an 
interferometric setup involving 2" optical paths [18j |. In 
this case the price to pay is the exponential growth of the 
number of optical paths and optical devices required for 
the implementation. Yet another way is based on the use 
of a single particle with 2" energy levels, where each level 
will embody a computational state of the register [(J . 



All these suggestions are inherently based on the rea- 
soning that usually a quantum algorithm consists only of 
transformations indiscriminately acting on all the qubits 
of a register (operators acting on all the states, for ex- 
ample H® n ) or specific state transformations (operators 
acting only on specific states, for example the phase- 
inversion of the searched state in the quantum search 
algorithm described in Sec. Ill A")) . However, the quantum 
average algorithm requires operators acting on specific 
qubits (for example the single-qubit phase-shifts). To 
give a clear picture of the differences between the two 
kinds of transformation above and operators acting on 
specific qubits, we consider the phase-shift stage in the 
quantum average algorithm. Two possible ways to simu- 
late it with classical processes can be distinguished. We 
can consider a serial sequence of transformations, each 
one equivalent to a single-qubit rotation (if we simulate 
the Hilbert space of n qubits in the above-mentioned 
ways, each transformation classically corresponds to an 
operator acting on 2™ _1 degrees of freedom of the classi- 
cal system 19]). However, we will lose the parallel com- 
putation characteristic, that is at the basis of this algo- 
rithm. The other way to classically simulate this stage 
is to consider a parallel application of all these transfor- 
mations. But it is easy to see that, considering again 
the ways mentioned above to simulate the Hilbert space 
of the quantum system, such a parallel application re- 
duces to a single transformation rotating one state of the 
classical system by an angle proportional to the average 
value (it also rotates the other states by different an- 
gles, but we are not interested in them). We thus need 
to know a priori the result of the algorithm. Therefore, 
the only non-trivial simulation changes radically the na- 
ture of the algorithm from parallel to serial computation. 
This highlights an intrinsic difference between the algo- 
rithms involving only global transformations or specific 
state transformations (polylocal class) and those which 
need transformations affecting only specific qubits (non- 
polylocal class). 
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III. NOISE EFFECTS ON THE 
REPRESENTATIVE ALGORITHMS 

Our model for noise is motivated by considerations typ- 
ical of static quantum chaos, where a register is assumed 
to be affected by individual, time-independent imperfec- 
tions on each qubit ■ Here we consider the possibility 
of an imperfect preparation of the state of the register by 
allowing each qubit to be in a mixed state. Intuitively, 
the loss of purity of the overall state can be expected to 
influence the behavior of the entanglement involved, if 
any, in a specific quantum algorithm. We therefore con- 
sider the initial state of the jth qubit (j = l,..,n) as 
given by the density matrix 



where Xj is the probability of finding the jth qubit in its 
ground state and 1 0) ^- is the ideal starting state. We name 
the source responsible for such an initial state as static 
noise. For our purposes, we do not need to identify the 
mechanism responsible for such imperfections to occur. 
This is a setup-dependent issue that will specialize our 
study. Nevertheless, we mention that if each Xj follows 
the Boltzmann distribution for a two-level system, this 
model for mixedness represents a qubit being thermally 
excited. Such an assumption is not at all unrealistic: the 
study of quantum algorithms in the presence of non-ideal 
preparation can be pragmatically relevant. For instance, 
solid-state implementations require the cooling of a regis- 
ter to very low temperatures, which may be experimen- 
tally demanding and quite unnecessary if a protocol is 
known to perform adequately with tolerable mixed ini- 
tial state. These considerations make our investigation of 
practical importance. We remark that this assumption 
of a static model for noise is only the first step toward a 
more complex and complete study of an important prob- 
lem in quantum information science. 



A. Polylocal class: Quantum search algorithm 

The simulation of noise effects can be archived by 
choosing a set of A/s and evaluating the probability of 
obtaining the searched state after each iteration of the 
Grover iterate. In order to fix the ideas and present our 
results in a clear way, we have considered a symmetrical 
case of all equal noise parameters Xj — A, Vj = 1, ...,n. 
This choice is reasonable for a spatially localized regis- 
ter, where all the qubits experience a noise mechanism of 
negligible strength fluctuations or are in touch with the 
same thermal environment. 

In Fig. [5] we present the results for a system of 2, 3 and 
4 qubits (panels (a) , (b) and (c) respectively) where we 
plot the probability P of finding a searched state against 
the number of iterations m of the Grover iterate and 
the noise parameter A. One can clearly see that the pe- 
riod of the oscillations between the two extremal states 



(a) 




FIG. 2: (Color online) Probability P to find a searched state 
(regardless of its specific instance) against the number of it- 
erations m and the noise parameter X in the quantum search 
algorithm acting on a noise-symmetric register of 2, 3 and 4 
qubits (panels (a), (b) and (c) respectively). 



involved in the algorithm (the equally weighted super- 
position |5) and the searched state) is the same for any 
value of noise. This property makes the protocol some- 
what robust against imperfections in the preparation of 
the register. In a situation where there is just a limited 
knowledge about the initial purity of the register, we do 
not have to make the protocol adaptive, as the period is 
left unchanged. This indirectly confirms the resilience of 
the timing of the protocol outlined in [2l| for different ap- 
proaches than ours. However, as the noise increases, the 
probability for the noisy algorithm to find the searched 
item is strongly affected. For any number of qubits larger 
than 2 and in absence of noise (i.e. for A = 1), P is 
not exactly 1 after the number of iterations correspond- 



Pi=Aj|0>j(0| + (l-A i )ll),-(l|, (4) 
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FIG. 3: (Color online) Normalized probability P norm to find 
the searched state (regardless of its specific instance) against 
the noise parameter A in the quantum search algorithm acting 
on a register of 2 (■), 3 (♦) and 4 qubits (k). 



ing to the first maximum (for 2 qubits it is well-known 
that P = 1 after just one iteration). We have therefore 
considered the normalized probability Pnorm = P/ Pideai , 
where Pideai is the probability to obtain the searched 
state when A = 1 and we have calculated P and Pideai 
after the number of iterations corresponding to the first 
maximum. The results are shown in Fig. [3] for the case 
of 2 (■), 3 (♦) and 4 qubits (A). The plot reveals that 
Pnorm (which is independent of the state that has to be 
found) scales as A" with n = 2, 3, 4, demonstrating a se- 
vere fragility of the scheme to static imperfections. This 
result can be easily generalized to an arbitrarily inhomo- 
geneous set {A.,} and to any dimension of the register. 
We find that 



-e- 



-(r^ H jj^Vv 



(5) 



FIG. 4: Circuit for the quantum average algorithm Q includ- 
ing the entangling gates for the preparation of the input state. 
Here, the states entering the circuit are given by the pj de- 
fined in Eq. Q. CNOT gates, controlled by the first qubit in 
the register, are shown. 



mance of the algorithm with Aj = 1 will not be precisely 
|s)(s|. However, the contributions from the orthogonal 
states will be negligible. Therefore, in the case of an 
inhomogenous set of A/s, results qualitatively analogous 
to those presented here should be expected. The fragility 
of the algorithm to this simple model for imperfections 
must be looked at in terms of the modification suffered 
by the entangled state "created" by the Grover iterate in 
the course of the protocol @ . 



B. Non-polylocal class: Quantum average 
algorithm 

In order to give a full-comprehensive analysis of the 
effects of noise in the quantum average algorithm, we 
explicitly include the steps required to create the entan- 
gled resource consumed during the computation. For this 
purpose, we consider the control-NOT (CNOT) gate Q 



with Pnorm the normalized probability of obtaining the 
searched state in the presence of generally asymmetric 
noise. This result can easily be understood by closely 
looking at the model for imperfections. Consider for in- 
stance the two-qubit case; the initial state of the reg- 
ister is written as p 12 = poo |00><00| + p i |01)(01| + 
Pro 1 10) ( 10 1 + pn 1 1 1 ) ( 1 1 1 . The probabilities pij are re- 
spectively poo = AiA 2 , poi = Ai(l- A 2 ), pio = (1 - Ai)A 2 
and pn = (1 — Ai)(l — A 2 ). Due to the linearity of quan- 
tum mechanics, we can study the evolution of each of the 
states involved in jOi 2 separately. The state |00)(00| is the 
initial state of the register in the ideal case of Ai. 2 = 1. 
The evolution of the system dictated by the algorithm 
is unitary so that the states present in the evolved state 
will remain orthogonal to each other. In particular, the 
algorithm will transform any other computational state 
into a state that will be orthogonal to the searched one 
resulting out of the evolution of |00)(00|. Therefore, the 
probability of obtaining the searched state after the right 
number of iterations is precisely poo = AiA 2 , as there will 
be no contribution from any other state. We can extend 
this proof to any number of qubits, therefore arriving to 
the result of Eq. |5]). Of course, for a number of qubits 
larger than 2, the ideal output state after the perfor- 



CNOT = 



'1 0\ 

10 

1 

v0 1 0/ 



(6) 



written in the two-qubit basis {|00> , |01> , |10> , |11», 
where the first qubit is the control and the second is 
the target of the gate. The GHZ-like state needed for 
the algorithm to work when Aj = 1 is obtained start- 
ing from the ideal initial state |0). The application of 
H to the first qubit only, followed by a set of CNOT's 
with the first qubit as the control and all the other qubits 
in the register as the targets, results in the required re- 
source. The complete protocol is sketched in Fig. 0] The 
preparation stage of the algorithm is addressed explicitly 
here because the static noise affecting the register has 
influences on the form of the entangled state used in the 
computation. By incorporating these preliminary steps, 
we provide a more complete understanding of these in- 
fluences. 

The choice of a meaningful figure of merit to compare 
the algorithm in the absence of noise with the noisy one is 
not straightforward. A number of difficulties arise during 
a careful analysis. One might naively consider the state 
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fidelity 

F = Tr(pidPnoise) = l(lp + \H PnoiseH\i) + )i (7) 

with pid — H\ip + )i(ip + \H the ideal final state of the first 
qubit just before its measurement and p no ise the corre- 
sponding mixed state in the presence of noise. Its ex- 
pression can be obtained analytically but, even in the 
simple case of symmetric noise, this is too lengthy to 
be reported here. State fidelity of the output states 
is often used as a significant parameter in the evalu- 
ation of the performances of a protocol. However, in 
what follows, we show that considering F as a figure 
of merit leads to the wrong conclusions. In order to 
fix the idea, for the calculations we have chosen the 
set {vi, 1*1,1/3} = {-0.775,0.25,0.675} with 6 = 0.0625. 
This allows us to give a clear picture of our results. Ob- 
viously, other choices are equally valid. 

The application of the algorithm in the presence of 
noise produces the plot shown in Fig. [5] For small values 
of A, i.e. when each qubit is prepared in a state close to 
|1), F is almost ideal (F = 0.95). This would lead us to 
conclude that the protocol is effective even for a prepara- 
tion orthogonal to the one designed for the algorithm to 
work. However, the conclusion is erroneous as revealed 
by immediately calculating the value |/u no j se |/0, with the 
pedex reminding us that the algorithm has been run in 
the presence of noise. By assuming A = we obtain 
\p-noise\/0 = 0.36 rather than the true value \p,\/0 = 0.80, 
showing a considerable discrepancy. On the other hand, 
at A = 0.1 we have F = 0.77 (lower then the one for 
A = 0) and an estimate |/x n |/# = 0.94, much closer to the 
actual value of Therefore, the state fidelity would 
lead us to erroneously privilege the first case over the 
second, which actually delivers a more faithful estimate 
of the average. 

A more significant performance parameter is given by 
the distance ratio D = (|/i no j se | ~ ImI)/^; which measures 
the distance between the true and estimated average in 
units of 9. Obviously, D ~ implies that the noise does 
not spoil the accuracy of the computation. By using D 
in the same situation considered above, we obtain much 
more faithful information, as shown in Fig. [6] (a), where 

F 
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FIG. 5: Fidelity F of the output state in a noisy quantum 
average algorithm against the noise parameter A, for v\ — 
-0.775, v 2 = 0.25, v 3 = 0.675 and = 0.0625. 



we can notice that, for A ~ 0.1, the algorithm gives a 
reliable evaluation of the ratio \p\/0- This might seem 
surprising at first sight. However, it is easy to recognize 
that this is simply a fortuitous case due to the depen- 
dence of the algorithm on the actual set of i/j's. Indeed, 
suppose we swap the value of v\ and v^: obviously the 
average value remains unchanged. However, we obtain a 
different plot for D (see Fig. [5] (b) ) . This is not a feature 
of the chosen figure of merit but an intrinsic character- 
istic of the quantum average algorithm. In a different 
way to the search protocol, here the qubits in the reg- 
ister play unequal roles. Indeed, in addition to carrying 
information about the first element of the set {i/j} after 
the phase-shift stage, the first qubit is also responsible for 
the information about the average value after the mea- 
surements on all the other qubits [22| . In practice, this 
may represent a problem: in a noise-asymmetrical set- 
ting, the noise affecting the last qubit to be measured 
is critical in determining the "quality" of the evaluated 
average. Ideally one would like to screen it from noise in 
order to have a more faithful estimate. 

A possible way to circumvent this problem is to con- 
sider a variation of the algorithm in which an enlarged 
register of N + 1 elements is used. The first qubit is mea- 
sured at the end of the algorithm while the rotations are 
performed on all the other N qubits. We assume that the 
first qubit is protected and in a pure state, while all the 
others are prepared in p/s. This situation is reminiscent 
of analogous investigations performed with respect to a 
different quantum algorithm [231 ]. Moreover, this scheme 
resembles the paradigm used in the model for determinis- 
tic quantum computation with one quantum bit [24| . Both 
the above cases showed that one pure state qubit singled 
out from a register prepared in a statistical mixture is 
sufficient to carry out several computational protocols. 
Our study reinforces such ideas and at the same time 
suggests an operative way to limit the effects of static 
imperfections. It is important to stress that by shield- 
ing the ruler qubit in such modified protocol, we want to 
effectively avoid the accumulation of noise effects rather 
than fix the mistakes occurred during the computation. 
This is different, in both motivations and strategy, from 



(a) (b) 

D D 




FIG. 6: (a): Distance ratio D of the output state of the 
quantum average algorithm against the noise parameter A, 
for the same set of values of Fig. [5] (b): Distance ratio D 
of the output state of the quantum average algorithm against 
the noise parameter A, for the same set of values of Fig. [5] 
but after v\ and v-2 have been swapped. 
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FIG. 7: Circuit for the modified quantum average algorithm 
with the introduction of the ruler qubit, which is prepared 
in a pure state that is not phase-shifted in the course of the 
protocol. All the remaining qubits enter the algorithm in 
Pi(j = 2,...,N + l). 



quantum error correction techniques [3, HE] . No redun- 
dant encodings is introduced, in our scheme, which are 
instead typical of error correction protocols. Each qubit 
in our modified scheme is a physical information carrier 
rather than a logical one encoded into the Hilbert space 
of a block of qubits. 

The logical circuit of the modified algorithm is shown 
in Fig. [3 where one can see that the required phase-shifts 
are now performed on all the qubits except the first in a 
register of N + 1 elements. The first qubit, also known as 
the ruler, has to be physically distinct with respect to all 
the others. This is in line with a scenario of distributed 
computation, where the quantum average algorithm was 
conceived: The register configuration can be that of a 
star graph with the ruler at the center and all the re- 
maining qubits occupying the outer vertices. Each qubit 
is connected to the ruler by classical and quantum chan- 
nels, needed to exchange the information acquired after 
the measurements and construct the entangled resource. 

Using the same set {vj} as before, we obtain the behav- 
ior of D shown in Fig. [H For A = we have l /x "°" e l = M 
so that for this value of A as well the algorithm works 
perfectly. Indeed, right before the application of the 
Hadamard and CNOT gates, the register is in the pure 



state |0> 1 



1 



N+l 



which then takes the form 



N+l 



V2 



(|0> 1 |1) 2 - ■ |1) 



N+l 



|l)i|0) 2 --|0) 



N+l/ 



(8) 



the 



This is still a GHZ-like state and thus of 
rect entanglement structure (formally N+1 
&xi n+i)- By performing the algorithm with 
we obtain the same probability cos 2 (^) (sin 2 (^)) that 
the ruler qubit is in |0) 1 ( 1 1) x ) . The dependence on the 
order of the values has disappeared, so that the use of the 
distance ratio in a modified protocol with a ruler qubit 
now allows us to perform a faithful assessment of a noisy 
algorithm. 

However, a further problem to address is the depen- 
dence of D on the average value. Here, we study the 
maximum amount of noise that the algorithm can toler- 
ate without affecting the ratio \(J,\/0, independently of the 
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FIG. 8: Distance ratio D of the modified quantum average 
algorithm with a ruler qubit against A and for the same set 
of values of Fig. [5] 



now be symmetrical with respect to A = i (correspond- 
ing to pj 's being completely mixed) , we decide to use the 
purity parameter t = |2A — 1 1 to quantify the strength of 
the imperfections. Obviously, for t = the noise will be 
maximum (all the input qubits will be completely mixed) 
while for r = 1 the input qubits will be in the pure state 
|0)^ or |1) ■ (in which case the ratio \n\/9 will be correctly 
evaluated, as we have shown before). The probability of 
obtaining |0) 1 in the presence of noise is evaluated to be 



Pn{t) = arccos[X^= (l 



,'N 
In, [ - 



>(-iy T 2i A 



(2i) 

(N—2i) 



with 




for even N 



for odd N 



(9) 



and 



perm 

. (10) 

Here, V is the permutation operator for the indices j 
and fc, which are in number of I and m respectively. For 
example 

A (i)- sin (y) sin (y) cos (y) + 



;vi\ . (v 2 \ . (v-i 

cos — sin — sin 



(11) 



It is straightforward to notice that Pjv(1) = cos 2 ^)- 
Rather than find the values of i/j/0 that maximizes 

|Pjv(t) - Pjv(1)| = llAWsel - | Mid 1 1 /0 f° r a set vame °f 
N, we can maximize the absolute value of the difference 
of the arguments in the inverse cosine as this is a mono- 
tonic function in the interesting range of values. We find 



that 



(-l)VM 



2i A (2i) 



(N — 2i) 



set v 3 (and therefore of the average value) and 6. As our maximum for z/j/0's all equal to a v. 
task is to check that \p\jB is 0(1), a precision of 0.5 will 
be considered as acceptable. As the behavior of D will 



_ l)*^(iV-2i)l 1S 

depending only 

on the number of qubits but otherwise independent of 



t. We have numerically calculated the value of u 
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FIG. 9: (Color online) Maximum value of \D\ against the 
purity parameter r for a register of 3 to 8 elements. In this 
plot, the number of qubits diminishes by one (starting from 
N = 8) in going from the top curve to the bottom one. 



for N — 3, ...,8 and report in Fig. [9] the corresponding 
\P N (r) - P N (1)\. 

The noise that the algorithm can tolerate while still 
giving a faithful estimate of \fi\/6 with a precision of 0.5, 
regardless of the values of Vj/6, is r ~ 0.90 and slightly 
depends on the number of qubits. The reason for this no- 
ticeable resilience to static imperfections is explained by 
analyzing the global density matrix of the register during 
the performance of the algorithm. We give an account of 
this robustness in the next section, which sheds light on 
the behavior of entanglement in the protocol itself. 



IV. ANALYSIS OF ENTANGLEMENT IN THE 
QUANTUM AVERAGE ALGORITHM 

The study of the entanglement behavior in the pres- 
ence of noise can help to understand the reasons for the 
resilience we have highlighted above. The fragility of the 
GHZ-like resource consumed in the algorithm may lead 
us to think that even simple static imperfections in the 
register would have dramatical effects on the quality of 
the computation. 

In analogy with what was done in the previous sec- 
tion, we consider the modified version of the quantum 
average algorithm with a ruler qubit. As the conditional 
phase-shifts used are local unitary operations, the degree 
of entanglement in the system is not modified by their 
application. Therefore, we can study the entanglement 
in the register before their action (of course the measure- 
ments performed at the end of the protocol will consume 
the entanglement). It is important to note that the en- 
tanglement dynamics are in contrast to those responsible 
for the speed-up in polylocal algorithms. In those cases, 
entanglement has to be created during the protocol and 
spread all over the register in order for the algorithms 
to outperform their classical analogs @, d, @. For the 
quantum average algorithm, a set of measurements pro- 
gressively breaks preconstituted quantum correlations. 

The focus of our interest is bipartite entanglement 
in the system. For this analysis we consider the 
[N + l)-qubit system as split into two subgroups. 



We reveal entanglement between bipartitions using the 
Peres-Horodecki negativity of partial transposition cri- 
terion 26]. Even though this test is not necessary and 
sufficient for revealing quantum correlations in the case 
of a general multipartite register, it is a useful tool for the 
present analysis. We first examine a three-qubit system 
and then generalize the results to any number of qubits. 

When dealing with three qubits, one can either con- 
sider the bipartite entanglement in the reduced state ob- 
tained by tracing out the degrees of freedom of one of the 
qubits (we call it the traced case) or look for the correla- 
tions between one qubit and the remaining two, consider- 
ing every possible permutation of qubit labels (we refer 
to this as the non-traced case). In the traced case, we 
have considered both the trace with respect to the ruler 
qubit (therefore studying the presence of bipartite entan- 
glement between two mixed qubits) and the trace with 
respect to one of the mixed qubits of the register (thus 
evaluating the quantum correlations of the remaining reg- 
ister qubit with the ruler). The result is that in both 
cases no bipartite entanglement is present. By studying 
the non-traced case, regardless of the configuration of the 
bipartitions, we have obtained that Vr ^ bipartite en- 
tanglement is present. These results are a reminder of 
the properties of GHZ states. General considerations for 
any number of qubits can be obtained by analyzing the 
static noise we have chosen. With this model, the den- 
sity matrix of the whole system, before the phase-shift 
stage, can be seen as an ensemble of density matrices of 
GHZ-like form. For instance, the global density matrix 
for the three-qubit system in the presence of noise is 



p = \ 1 \*\GHZ$){GHZ& 
+ X^l- \ 2 )\GHZ$)(GHZ 01 
+ (l-X 1 )X 2 \GHZg ) )(GHZg ) \ 
+ (l-\ 1 )(l-\ 2 )\GHZ[f)(GHZff\ 



r(2)| 



(2) 



(12) 



where we have introduced the set of generalized GHZ-like 
states 

|G#4^ -1(|000) + |111)), 



\GHZ, 



(2) 



1 



01 y/2 



\GHZ 



(2) 



GHZ[f) = -==(|011) + 1 100)). 



V2 
1 



(|001) + |110)), 
(|010) + |101)), 



(13) 



In general, for an iV-qubit system, the global density ma- 



trix will be the sum of projectors \GHZ ( ^\){GHZ ( ^ ) 



with 



\GHZ\% } ) = -J=(|0 ai -a N ) + |1(1 - oi)-(l - on))) 

(14) 

and {ai} the ordered sequence of digits of a binary 
number between and 2^—1. The coefficent of each 
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\GHZg> y ) is = n,\ ? '(l " A,) 1 " 01 - When deal- 

ing with an initial state |Gi/z|^), the protocol gives 
an estimate of the absolute value of a modified average 
jEt = jf J2f =1 (— \) aj Vj. In presence of a not-too-severe 

noise, the dominating coefficents ci N \ are those with 
fewer l's in the set {ai}. The corresponding estimate 
of the average will be close to the actual one. On the 
other hand, the errors that result from the states with 
a large estimate discrepancy (those with a number of 
l's in the set {a{\ close to N/2) will be damped by the 
corresponding small coefficents cj^j. This explains the 
reason why the algorithm is robust. This same reasoning 
can be applied to explain the entanglement behavior of 
the register during the performance of the algorithm. In- 
deed, the state of the register can be seen as an ensemble 
of GHZ-like states that will maintain the corresponding 
characteristics even for large values of A. 

We have also analyzed the algorithm under the effects 
of a different model of static noise. Assuming that the 
initial state of the register is 

p = f \GHZi$. )(GHZi£>. \ + ^ I (15) 

(i.e. we are now considering white noise) the algorithm 
becomes fragile also for f ~ 1 as the state of the register 
can no longer be seen as an ensemble of GHZ-like states. 
This is further proof that the GHZ-like nature of entan- 



glement is a pre-requisite for the efficiency of this specific 
algorithm. 



V. REMARKS 

We have studied a representative of the non-polylocal 
class of algorithms (in contrast with the polylocal class 
considered by Meyer Q). We have found that a model 
of static noise preserving the GHZ-like nature of the en- 
tanglement in the register has no dramatic effects on the 
algorithm. On the other hand, if this specific entangle- 
ment structure is lost, the efficiency of the algorithm in 
term of accuracy is compromised. These results imply 
that the GHZ-like nature of the entanglement has a fun- 
damental role in this specific protocol. The dependence 
on other classes of entanglement in other non-polylocal 
algorithms deserves further investigation and will be the 
subject of future study. 
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